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Abstract 

We study invariants under gauge transformations of linear partial differential op- 
erators on two variables. Using results of BK-factorization, we construct hierarchy of 
general invariants for operators of an arbitrary order. Properties of general invariants 
are studied and some examples are presented. We also show that classical Laplace 
invariants correspond to some particular cases of general invariants. 
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1 Introduction 

Laplace invariants a = c — ab — a x and b = c — ab — b y , introduced for bivariate 
hyperbolic operator of the second order 

d x d y + ad x + bdy + c, (1) 

are quite well-known and have been studied by many researchers (see, for instance, 
[2j and bibl. herein). Their importance is due to the classical theorem: 

Theorem 1 Two operators of the form are equivalent under gauge transformations 
iff when their Laplace invariants coincide pairwise. 

These invariants are used extensively in the theory of integr ability ([Hj, [2]) |1 lj ■ 
|9|.[T0] and others). Obviously, Laplace invariants can be regarded as factorization's 
"remainders" for the initial operator Q: 



d x d y + ad x + bd v + c 



(d x + b){dy + a) - ab- a x + c , . 

(d y + a)(d x + b)-ab-b y + c U 

and this presentation is a first step of Darboux-Laplace transformations used for 
construction of solution to non-factorizable linear partial differential equations in two 
independent variables. 

To generalize construction of invariants to the case of an arbitrary order operator, 
a factorization algorithm for such an operator is needed. Recently two factorization 
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algorithms for arbitrary order bivariate LPDOs have been published. In p], the fac- 
torization algorithm called Hensel descent is presented, which is a modification of well- 
known Hensel lifting algorithm (see, for instance, [H) and which allows to factorize 
operators from the ring D = Q(x,y)[d x ,d y ]. In pQ, absolute factorization is presented 
for operators whose coefficients are arbitrary smooth functions of two variables. The 
idea to connect the results produced by this method (called now BK-factorization) and 
Darboux-Laplace transformations has been briefly discussed in 5 . In present paper we 
construct general invariants using BK-factorization procedure which can be described 
as follows. Factorization problem is regarded in the form 

A n = a jk dld k y = ( Pl d x +p 2 d y + P3 )( Yl I'j^i) (3) 

j+k<n j+k<n 

with assumption that a n Q ^ (locally) and can be taken as o n o = 1 without loss of 
generality. It is shown that for any simple root u of characteristic polynomial 

n 
k=0 

with n > 2, there exists unique factorization of the form (|3J): 

n-2 

A = LB + C, L = d x -ujd y +p, C = ^2cjdi. (4) 

and the conditions of factorization are written out explicitly. In Section 2 we use 
these results for construction of general invariants of an arbitrary LPDO and study 
some of their properties while in Section 3 the whole hierarchy of general invariants is 
described. In Section 4 a few instructive examples are given. In particular, it is shown 
that both Laplace invariants can be constructed as a simple particular case of general 
invariants. In Section 5 brief discussion of the results obtained is presented. 



2 General invariants and semi- invariants 

The conditions of factorization produced by BK-factorization have following form: 
for second order operator 

A 2 = a jkdi d y = (Pidx + P2d y + P3)(P4d x + p^dy + p 6 ) 

j+k<2 

the only factorization condition is 



aoo = L(p G ) +P3P6, (5) 

for third order operator 
A 3 = ^ a jkd J x dy = (pid x + p 2 d y + pz){pidl + p 5 d x d y + p 6 dy + p 7 d x + p 8 d y + p 9 ) 

j+k<3 

there are two factorization conditions: 

|a i = L{ps) +P3P8+P2P9 ^ 

\a Q0 = L(p 9 ) +p 3 P9 
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and so on for arbitrary n. Now we can give explicit representation of a LPDO whose 
characteristic polynomial has at least one simple root, in the form of factorization with 
remainders: 



A 2 = (pid x + p 2 d y + P3){p4d x + p 5 d y + p 6 ) - l 2 
with remainder l 2 = aoo — L(pq) + P3Pe] 

A 3 = (pid x + p 2 d y + P3){pa91 + p 5 d x d y + p 6 d y + p 7 d x + p 8 d y + p 9 ) - l 3 d y - Z 31 , 

with two remainders Z 3 = a i - (pid x + p 2 d y +P3)Ps -P2P9, and I31 = a 00 - (pid x + 
P%9y + P3)P9; and so on for arbitrary order of LPDO. 

Definition The operators A, A are called equivalent if there is a gauge transforma- 
tion that takes one to the other: 

Ag = e^A(e^g) = A v g. 

Note that 

(d x )<p = dx + <Px, {9y)<p = 9y + tfy 

and a gauge transformation is an algebra automorphism. Therefore A and have 
the same characteristic polynomial and the factorization (jlj carries over: 

Afp = L^B^ -\- C^. 

It follows that both the characteristic polynomial and the leading nonzero coefficient 
of the remainder term are invariants. Number of remainders varies for operators of 
different orders, and also their properties are different. In order to demonstrate it, let 
us formulate the following 

Lemma 1 For an operator of order 2, its remainder l 2 is its invariant under the 
equivalence transformation, i.e. 

h = h- 

For an operator of order 3, its remainder I3 is its invariant, i.e. 

^3 = h, 

while remainder I31 changes its form as follows: 

^31 = ^31 + hPy 

► Introducing notations A 2a = A 2p - l 2 , A 3a = A 3p - l 3 d y - I31, f = e v , we 

get 

A 2a = r x A 2a o f = r\A 2p -i 2 )of = A 2p - r\h) <> / = A 2p - 1 2 . 

and 

-43a = f~ 1 A3 a of = f- 1 (A 3p -l 3 d y -l 3 i)of = A 3p - f' 1 (hi+l 3 d y )of = A 3 p-l 3 dy-l 3 <p y -l 31 M 
Corollary 1: If Z3 = 0, then Z31 becomes invariant. 

That is the reason why we call l 3 \ further semi-invariant. 



3 



Corollary 2: If h 7^ 0, it is always possible to choose some function / : hi = 

h<Py + ^31 = 0. 

Note that for second order operator, if its invariant l 2 = then operator is fac- 
torizable while for third order operator two its invariants have to be equal to zero, 
^3 = hi = 0. On the other hand, if operator of third order is not factorizable we can 
always regard it as an operator with only one non-zero invariant. Of course, all this is 
true for each simple root of characteristic polynomial, so that one expression, say, for £3, 
will generate three invariants in case of three simple roots of corresponding polynomial. 



3 Hierarchy of invariants 

As it was shown above, every general invariant is a function of a simple root to of the 
characteristic polynomial and each simple root provides one invariant. It means that 
for operator of order n we can get no more than n different invariants. Recollecting 
that BK-factorization in this case gives us one first order operator and one operator 
of order n — 1, let us put now following question: are general invariants of operator of 
order n — 1 also invariants of corresponding operator of order n? 

Let regard, for instance, operator of order 3 
A 3a = AiA 2a - hd y - hi = Ai(A 2p - h) - hd y - hi = A x A 2p - l 2 A x - l 3 d y - hi 
and 

A?,a = A\A 2v — l 2 Ai — hdy — hi = AiA 2p — l 2 Ai — hd y — hi, 

i.e. l 2 is also invariant of operator A 3a . Let us notice that general invariant h = h(w^) 
is a function of a simple root u/ 3 ) of the polynomial 

Vz{z) = a 30 z 3 + a 2 iz 2 + ai 2 z + a 3 

while general invariant l 2 = l 2 (u^) is a function of a simple root u/ 2 ) of the polynomial 

^2 (z) = Paz 2 +P5Z+P6 

with P4, j>5, pq given explicitly for co = a/ 3 ). In case of all simple roots of both 
polynomials Vz(z) and lZ 2 (z), one will get maximal number of invariants, namely 
6 general invariants. Repeating the procedure for an operator of order n, we get 
maximally n\ general invariants. In this way for operator of arbitrary order n we can 
construct the hierarchy of its general invariants 

l>ni ln—li ■■■jh 

and their explicit form is given by BK-factorization. As to semi-invariants, notice that 
an operator of arbitrary order n can always be rewritten in the form of factorization 
with remainder of the form 

Indx + in,l^ _1 + - + L,k-1, k<n 

and exact expressions for all U are provided by BK-factorization procedure. The same 
reasoning as above will show immediately that l n is always general invariant, and each 
ln,k-i i s io-th semi- invariant, i.e. it becomes invariant in case if l n ,k-i = 0, Vi < iq. 
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Example 1 Let us regard a third order hyperbolic operator in the form 

C = 030^ + a,2\d 2 j.d y + a\ 2 d x dy + ao3<9y 3 + terms of lower order (7) 

with constant high order coefficients, i.e. a^- = const V % + j = 3 and all roots of 
characteristic polynomial 

a 30 w 3 + a 21 u 2 + a 12 u + a 03 = V 3 (uS) 

are simple and real. Then we can construct three simple independent general invariants 
in following way. Notice first that in this case high terms of ((7J) can be written in the 
form 

(aid x + Pid y )(a 2 d x + (3 2 d y )(a 3 d x + (3 3 d y ) 

for all non-proportional otj, (3i and after appropriate change of variables this expression 
can easily be reduced to d x d y (d x + d y ). Let us introduce notations d\ = d x , d 2 = 
dy, 9% = d\ + d 2 = dt, then all terms of the third and second order can be written out 
as 

Cijk = (di + a,i)(dj + aj)(d k + a k ) = didjd k + a k didj + ajdid k + aidjd k + 
+ (dj + aj){a k )di + (<% + a,i)(a k )dj + (di + ai)(aj)d k + (di + a i )(d j + aj)(a k ) 

with 

020 = a 2 , a 02 = oi, an = ai + a 2 + a 3 

and Cjjfc = C — C^ k is an operator of the first order which can be written out ex- 
plicitly. As it was shown above, coefficients of Cij k in front of first derivatives are 
(semi) invariants and therefore, any linear combination of invariants is an invariant it- 
self. Direct calculation gives us three simplest general invariants of the initial operator 
C: 

hi = a 2jX — ai jV , Z32 = 03,4/ — a 2,t, hi = 03^ — a\ yt . 

Lemma 2 General invariants l 2 \, l 32 , hi are ah equal to zero iff operator C is 
equivalent to an operator 

l = d x d 2 d z + hdt + b 2 d 2 + c , (8) 

i.e. 3 function / : f~ l C o f = L. 
► Obviously 

f-\d x dyd t ) 0f=(d x + (lOg f) x )(dy + (lOg /)„)($ + (tog/) t ) 

for any smooth function /. Notice that it is the form of an operator C^ k and introduce 
a function / such that 

ai = (log/) a ., a 2 = (logf) y , a 3 = (log/) t . 
This system of equations on / is over-determined and it has solution /o iff 

a2,x ~ ai, y = 0, a 3 , y - a 2 ,t = 0, a 3)X - a^ t = 0, (9) 

i.e. hi = l 32 = hi = 0. 

Indeed, if C is equivalent to ©, then 020 = 002 = an = and obviously 
hi = I32 = hi = 0. ■ 

Note that the statement of the Lemma 2 is weaker than that of Theorem 1. 
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Example 2 First two equations of © determine solution uniquely, up to a constant, 
and therefore determine 

<%/ = ipi + d 2 )f = ai + a 2 . 

However the compatibility conditions are satisfied if d 2 a\ = d\a 2 and (Z3 = a\ + a 2 + c, 
c constant. In general, the compatibility conditions only tell that two equations of 
are satisfied, but the third may be off by a constant. For instance, two constant 
coefficient operators 

d x dy(d x + d y ), d x d y (d x + d y + 1) 
have the same general invariants (all zero) but are not equivalent. 

4 Further examples 

Example 3 Let us regard hyperbolic operator in the form 

d xx - d yy + a w d x + a 01 d y + a 00 , (10) 
i.e. a2o = 1, an = 0, a$ 2 = —1 and to = ±1, C = d x — ud y . Then l 2 takes form 
h - a o - £( ^ Ql ° ~ aQ1 ) _ ^ Q io ~ Q oi ^"io + a oi 

which yields, for instance for the root u> = 1, to 

7 - n n a ^ ~ a oi \ a io ~ Q oi _ „ r a aip - a i Q io ~ a oi 
t2 - aoo - 5 1 4 ~ a °° ~ ( dx ~ d v){ 2 ' 4 

and after obvious change of variables in (|10|) we get finally first Laplace invariant a 

l 2 = c — d x a — ab = a, 

where 

«io — 001 , 010 + a oi 
a= , 0= , c = a 00 . 

Choice of the second root, u = —1, gives us the second Laplace invariant b, i.e. Laplace 
invariants are particular cases of the general invariant so that each Laplace invariant 
corresponds to a special choice of u. 

Example 4 Let us proceed analogously with an elliptic operator 

d xx + d yy + a w d x + a 01 d y + a 00 , (11) 
then lo = ±i, £ = d x — ujd y and 

h = aoo + {d x T ^ y ){ ) + 1 

where choice of upper signs corresponds to the choice of the root u = i and choice of 
lower signs corresponds to w = —i. 
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Example 5 Now let us regard a third order operator 

B = d 2 d y + d x dy + a n d x d y + a 10 d x + a m d y + a o, (12) 

with CI30 = ao3 = a2o = ^02 = 0, «2i = 012 = 1 and all other coefficients are 
functions of x, y. Then (semi) invariants 

h = x a u - o i and Z 31 = c^aio - a o 

have very simple forms and gives us immediately a lot of information about the prop- 
erties of operators of the form (|12jl. for instance, these operators are factorizable, i.e. 
has zero invariants Z3 = Z31 = 0, iff 

a u = J a 01 dx + fi(y), a w = J a m dx + f 2 (y) 

with two arbitrary functions on y, fi(y) and f2{y)- 



5 Summary 

Now, that construction of Laplace invariants has been generalized to the case of arbi- 
trary order and arbitrary type of LPDO, the next obvious step would be to generalize 
of Darboux-Laplace transformations to this class class of operators, to construct gen- 
eralized Todda lattice, etc. It would require a deeper study of intrinsic properties of 
general invariants. For instance, beginning with operator of order 4, their maximal 
number is bigger then number of coefficients of a given operator, 

(n + l)(n + 2) , w 

^ ^ '- <n\ Vn > 3. 

It means that general invariants are dependent on each other and it will be a chal- 
lenging task to extract the subset of independent general invariants, i.e. basis in the 
finite space of general invariants. Another important step would be to extract a subset 
providing full analog of Theorem 1 for arbitrary order operators. 

Already in the case of three variables, the factorization problem of a corresponding 
operator and also constructing of its invariants becomes more complicated, even for 
constant coefficients. The reason of it is that in bivariate case we needed just to 
factorize leading term polynomial which is always possible over C. It is not the case 
for more then 2 independent variables where a counter-example is easily to find: 

Example 6 A polynomial 

x 3 + y 3 + z 3 — 3xyz = (x + y + z)(x 2 + y 2 + z 2 — xy — xz — zy) 

is factorizable but it is easy to prove that x 3 + y 3 + z 3 is not divisible by a linear 
polynomial 

ax + 0y + 72; + 5 

for any complex coefficients a, (3, 7, S. Thus, some non-trivial conditions are to be 
found for factorization of polynomials in more then two variables. 
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